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Measurement of magnetic fields on the few-hundred-kilometer length scale is significant for a variety of geo- 
physical applications including mapping of crustal magnetism and ocean-circulation measurements, yet avail- 
able techniques for such measurements are very expensive or of limited accuracy. We propose a scheme for 
remote detection of magnetic fields using the naturally occurring atomic-sodium-rich layer in the mesosphere 
and existing high-power lasers developed for laser guide-star applications. The proposed scheme offers dramatic 
reduction in cost, opening the way to large-scale magnetic mapping missions. 



I. INTRODUCTION 



Measurements of geomagnetic fields are an important tool 
for peering into the earth's interior, with measurements at 
differing spatial scales giving information about sources at 
corresponding depths. Measurements of fields on the few- 
meter scale can locate buried ferromagnetic objects (e.g. un- 
exploded ordnance or abandoned vessels containing toxic 
waste), while maps of magnetic fields on the kilometer scale 
are used to locate geological formations promising for mineral 
or oil extraction. On the largest scale, the earth's dipole field 
gives information about the geodynamo at depths of several 
thousand kilometers. Magnetic measurements at intermediate 
length scales, in the range of several tens to several hundreds 
of kilometers likewise offer a window into important scientific 
phenomena, including the behavior of the outer mantle, the 
solar-quiet dynamo in the ionosphere UJ, and ionic currents 
as probes of ocean circulationlJ], a major actor in models of 
climate change. 

To avoid contamination from local perturbations, magnetic- 
field measurements on this length scale must typically be 
made at a significant height above the earth's surface. Though 
magnetic mapping at high altitude has been realized with 
satellite-borne magnetic sensors l[3l0]|5i, the great expense of 
multi-satellite missions places significant limitations on their 
deployment and use. Here we introduce a high-sensitivity 
method of measuring magnetic fields with 100-km spatial res- 
olution without the cost of spaceborne apparatus, by exploit- 
ing the naturally occurring atomic sodium layer in the meso- 
sphere and the significant technological infrastructure devel- 
oped for astronomical laser guide stars (LGS). This method 
promises to enable creation of global sensor arrays for contin- 
uous mapping and monitoring of geomagnetic fields without 
interference from ground-based sources. 
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FIG. 1 : Fluorescent detection of magneto-optical resonance of meso- 
spheric sodium. (Diagram not to scale). Circularly polarized laser 
light at 589 nm, modulated near the Larmor frequency, pumps atoms 
in the mesosphere. The resulting spin polarization (pictured as in- 
stantaneously oriented along the laser beam propagation direction) 
precesses around the local magnetic field. Fluorescence collected by 
a detection telescope exhibits a resonant dependence on the modula- 
tion frequency. 



II. PRINCIPLE OF THE METHOD 

The measurement we envisage is closely related to the tech- 
niques of atomic magnetometry, appropriately adapted to the 
conditions of the mesosphere. The principle of the method 
is to measure the magnetic field using spin precession of 
sodium atoms by creating atomic spin polarization, allowing 
it to evolve coherently in the magnetic field, and determin- 
ing the post-evolution spin state. Conceptually, these pro- 
cesses are distinct and sequential, although in practice they 
may occur simultaneously and at different times for differ- 
ent atoms. Preparation of spin-polarized mesospheric sodium 
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atoms is achieved by optical pumping, as was proposed for 
mesospheric sodium in the seminal paper on sodium LGS by 
Happer et al. 1 6 1 . In the simplest realization, the pumping laser 
beam is circularly polarized and is launched from a telescope 
at an angle approximately perpendicular to the local magnetic 
field, as shown in Fig. [T] In the presence of the magnetic 
field, the transverse polarization generated by optical pump- 
ing precesses around the magnetic-field direction at the Lar- 
mor frequency. In order to avoid "smearing" of the atomic 
polarization due to this precession, the optical-pumping rate 
is modulated near the Larmor frequency, as first demonstrated 
by Bell and Bloom 17|. When the modulation frequency and 
the Larmor frequency coincide, a resonance results, and a sub- 
stantial degree of atomic polarization is obtained. The atomic 
polarization in turn modifies the fluorescence from the sodium 
atoms, which can be detected by a ground-based telescope. 
This allows the sodium atoms to serve as a remote sensor of 
the magnitude of the magnetic field in the mesosphere, i.e. as 
a scalar magnetometer. By stationing lasers and detectors on a 
few-hundred-kilometer grid, a simultaneous map of magnetic 
fields may be obtained; alternatively, the laser and detector 
can be mounted on a relocatable stable platform such as a ship 
or truck to facilitate magnetic surveying. In contrast to both 
ground-based and satellite-based measurements, the platform 
is not required to be magnetically clean or quiet. 

The magneto-optical resonance manifests itself as a sharp 
increase in the returned fluorescence for the D2 line of 
sodium, or a decrease for the Dl line, as a function of the 
modulation frequency. The sodium fluorescence may be cap- 
tured either by a single-channel photodetector or by an imag- 
ing array (e.g., a CCD camera); the latter allows additional 
background discrimination and the possibility of separately 
analyzing fluorescence signals from different altitudes if the 
imaging and laser-beam axes are not parallel. 



III. MEASUREMENT PARAMETERS 

The magnetometric sensitivity of this technique is governed 
by the number of atoms involved in the measurement, the co- 
herence time of the atomic spins, and the fraction of the to- 
tal fluorescence intercepted by the detector, each of which we 
seek to maximize. A detailed discussion of mesospheric prop- 
erties relevant to LGS has recently been given [8|. Briefly, 
the altitude of the sodium layer is ~90km, its thickness is 
~10km, and sodium is present at a typical temperature of 
180 K and a number density of around 3 x 10^ m^'^. This 
density is low by vapor-cell standards, but the interaction vol- 
ume and hence the atom number can be quite large, limited 
chiefly by available laser power The coherence time is lim- 
ited primarily by collisions with other atmospheric molecules 
and secondarily by atom loss from the region being probed 
(e.g., due to diffusion or wind). A velocity-changing collision 
occurring after an atom is pumped typically removes the atom 
from the subset of velocity classes which are near-resonant 
with the laser light; as a consequence, these collisions result 
in an effective decay of spin polarization. Moreover, spin- 
exchange collisions of sodium atoms with unpolarized para- 



magnetic species in the mesosphere, predominantly O2, result 
in a randomization of the electron spin, and therefore also lead 
to decay of sodium polarization. To our knowledge, the spin- 
exchange cross-section of oxygen with sodium has not been 
measured, and in fact it is anticipated that measurement of 
mesospheric coherence times will more tightly constrain this 
important LGS parameter However, its magnitude can be es- 
timated from other known spin-exchange cross-sections, lead- 
ing to an expected spin-damping time on the order of 500 /is. 
The fraction of intercepted fluorescence is determined by the 
solid angle subtended by the detection telescope and the an- 
gular emission pattern of the fluorescing atoms; for a 1-m^ 
telescope and isotropic emission, the fraction is approximately 
10^^^ when the detector is directly below the fluorescing sam- 
ple. 

Sodium atoms in the mesosphere are distributed both in po- 
sition and in velocity. They can be conceptually divided into 
velocity classes and spatial volumes, in each of which the col- 
lection of atoms functions approximately as an independent 
magnetometer. The width in velocity Av of an individually 
addressable velocity class is set by the Doppler effect and the 
natural linewidth 70 « 27r x 10 MHz of the sodium transi- 
tion. We can thus approximately calculate the sensitivity by 
calculating the evolution of a single sub-volume and veloc- 
ity class; we then determine the total fluorescence by sum- 
ming over sub-volumes and velocity classes. This approxi- 
mate treatment is valid for length scales larger than the typical 
distance traveled by atoms in one coherence time, and for ve- 
locity classes that are approximately uncoupled. The latter 
condition requires that the laser intensity be low, so that radia- 
tion pressure is small, and that there be negligible probability 
of an atom begin transfered from one near-resonant velocity 
class to another by velocity-changing collisions. 

The atom number can be maximized either by increasing 
the interaction volume or by including more velocity classes. 
Since the most probable atomic velocity component along the 
laser beam propagation direction is zero, broadening the laser 
spectrum to include velocity classes away from zero velocity 
will offer diminishing returns as spectral widths approaching 
the Doppler linewidth (~1 GHz) are reached. By contrast, de- 
focusing the laser beam to illuminate a larger region in space 
suffers no limitation other than the amount of laser power 
that can be supplied at a reasonable cost. For this reason, as 
well as for simplicity of implementation, defocusing the laser 
beam appears advantageous. If the observed relaxation due to 
velocity-changing collisions is significantly more rapid than 
that due to spin-exchange collisions, however, it may prove 
beneficial to broaden the laser spectrum to a width comparable 
to the Doppler linewidth, so that velocity-changing collisions 
will no longer cause loss of polarization but will merely cause 
atoms to be pumped and probed in different velocity classes. 



IV. RESULTS 

One expects on intuitive grounds that the optimum laser in- 
tensity resonant with a single velocity class should be such 
that the characteristic rate of optical pumping Tp = 7o//2/sat 
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FIG. 2: Calculated magneto-optical resonance profiles for meso- 
spheric sodium.The resonances shown correspond to the D2 F=2 
F'=3 (upper curve, blue diamonds) and Dl F=2 F'=l (lower 
curve, green circles) sodium lines. The Dl curve has been multiplied 
by a factor of three to improve visibility. Symbols are the results 
of numerical calculations, and solid lines are Lorentzian fits to these 
results. Calculations are for an intensity I — 1 W/m^ and a modu- 
lation duty cycle of 20% with a detector collection area of 1 m^. 



(where / is the laser intensity and /gat ~ GOW/m^ is the 
saturation intensity of the sodium cycling transition) is on the 
same order as the decay rate of atomic polarization. For higher 
intensities, optical pumping will "reset" the precession before 
its coherence time has been fully exploited, while use of lower 
intensities sacrifices signal without improving the coherence 
time. We have performed a detailed ground-state density- 
matrix analysis of spin precession and optical pumping on the 
Dl and D2 transitions of sodium. In this analysis, we em- 
ploy a circularly polarized pump laser beam oriented at right 
angles to a magnetic field of 0.5 G, with a spin-exchange col- 
lision time of 500 /is and a velocity-changing collision time 
of 200 /is. Since we employ optical intensities substantially 
lower than the saturation intensity, we expect the ground-state 
method to be accurate; we have also performed calculations of 
the resonance contrast for selected parameters using the full 
(ground and excited-state) optical Bloch equations. For a va- 
riety of settings of the pump light intensity and of the duty 
cycle (defined as the duration of a single pumping light pulse 
divided by the period of the modulation), we calculate reso- 
nance spectra as functions of the laser modulation frequency. 
Sample spectra for the laser tuned to the D1F = 2^F' = 1 
and the D2 F — 2 F' = 3 lines are shown in Fig.|2] From 
the width and peak height of these spectra, as well as the op- 
tical shot noise of the detected fluorescence, we calculate the 
magnetometric sensitivity. We assume that the noise is dom- 
inated by the fundamental shot noise and therefore neglect 
technical noise due to the photometric measurement. Con- 
tour plots of the sensitivity are shown as functions of the duty 
cycle and the laser intensity in Fig.[3]for the Dl and D2 lines. 

In these calculations, a fixed number of velocity classes is 
considered (specifically fifteen, or five per excited-state hy- 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Duty Cycle 



FIG. 3: Contour plot of calculated magnetometer sensitivity as a 
function of pumping duty cycle and intensity.Sensitivity is calcu- 
lated including the hyperfine effect for the D2 F=2 F'=3 and Dl 
F=2^ F'=l lines, using a detector area of 1 rn^, a spin-exchange 
collision time of 500 fis, and a velocity-changing collision time of 
200 fj,s. Contours are logarithmically spaced at intervals of one oc- 
tave. 



perfine component), and the effective laser beam size is ad- 
justed to maintain a constant launched laser power of 20 W 
(typical of the latest generation of LGS lasers), which results 
in different numbers of participating atoms for different laser 
intensities. The optimum sensitivity of 1.2nT/\/Hz occurs 
on the Dl transition, at a duty cycle of 25% and pump inten- 
sity during each pulse of 1.0 W/m^. The optimum intensity 
coiTesponds to an effective laser beam diameter in the meso- 
sphere of around 5 m. This sensitivity is approximately one 
order of magnitude worse than the limit set by quantum spin- 
projection noise for the sodium atoms, presumably as a result 
of hyperfine structure and Doppler broadening. The optimum 
on the Dl line offers superior magnetometric sensitivity in 
part because the Dl resonances are dark, i.e., they result in 
a reduction of fluorescence, so that the photon shot noise is 
smaller and broadening of the resonance is reduced. If techni- 
cal rather than fundamental noise sources dominate, then the 
D2 resonance may be preferable for its larger signal size. 

We note also that the time scale of velocity-changing colli- 
sions is comparable to the relaxation time of the atomic spins, 
so that it is likely for a given atom to interact more than 
once with the laser light while remaining in a single velocity 
class. Thus the excited-state hyperfine structure is resolved 
throughout the pumping and probing process, and proximity 
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to a given ground-to-excited hyperfine transition can strongly 
influence the light-atom interaction[9J. For the same rea- 
son, magneto-optical resonances involving higher polariza- 
tion multipoles such as alignment (which can be prepared by 
pumping with linearly polarized light) are expected to be ob- 
servable in mesospheric sodium. 

V. ADDITIONAL EFFECTS 

The magneto-optical resonance linewidth in mesospheric 
sodium is broader than in typical vapor-cell magnetome- 
ters; as a result, several effects important in vapor cells are 
less significant for mesospheric-sodium measurements. The 
quadratic Zeeman shift, for instance, leads to a splitting of 
the resonance into multiple resonances spaced by ~150Hz. 
However, since the width of the resonance at the optimum 
is around 5 kHz, this splitting will merely result in a small 
(and calculable) distortion of the lineshape. The natural inho- 
mogeneity of the geomagnetic field will also affect the mea- 
surement, since the resonant frequency varies with altitude. 
Variation of the earth's field over the 10-km height of the 
sodium layer is on the order lfTOl of 100 nT, corresponding to 
about 700 Hz in Larmor frequency. Thus, although the mod- 
ulation pumping laser cannot perfectly match the resonance 
frequency throughout the height of the mesosphere, this inho- 
mogeneous broadening is again small compared to the width 
of the resonance at the sensitivity optimum. Consequently, 
we expect the effects of both the natural magnetic gradient 
and the quadratic Zeeman shift to be small. Temporal varia- 
tions of the magnetic field in the mesosphere are, in principle, 
merely part of the signal being measured, and not an instru- 
mental limitation. However, large enough fluctuations could 
make it difficult to track the resonance frequency. We take as 
a likely upper bound for the magnetic fluctuations on time- 
scales of 1 s to 100 s the typical observed value at the earth's 
surface under ordinary conditions of around 1 nT. As this is 
again substantially smaller than the resonance linewidth, we 
expect that except during magnetic storms, it should not be 
difficult to keep the laser modulation frequency on resonance. 
Variations of the height and density of the sodium layer itself 
are an additional practical concern. A realistic measurement 
will require reducing sensitivity to such variations through 
comparison of on-resonant and off -resonant signals, either by 
temporally dithering the modulation frequency or by employ- 
ing spatially separated pump beams with different modulation 
frequencies. 

A further deviation of the real experiment from the ideal- 
ization embodied in the calculations comes from turbulence 
in the lower atmosphere, which causes random phase shifts 



in distinct transverse patches of the laser beam. In the typi- 
cal LGS application, the far-field diffraction (or speckle) pat- 
tern in the mesosphere from these low-altitude phase patches 
consists of elongated filaments whose individual lateral size 
is set by the numerical aperture of the laser launch telescope, 
but whose collective extent is governed by the size of the at- 
mospheric patches|6|. Fluctuation of these filaments in time 
results in undesirable effects including random changes in 
pump-laser intensity in the mesosphere, motion of the illu- 
minated column, and variation of the returned fluorescence. 
The relatively low intensity and large beam area indicated 
by our magnetometry calculations make such filamentation 
a lesser concern, however, since the beam diameter may be 
kept within the Fried length ifTTIl of around 0.1m in the most 
turbulent region of the lower atmosphere. Although lensing 
and beam-steering due to atmospheric variations will prevent 
precise fine-tuning of laser intensity, we do not anticipate that 
such variations will strongly affect the sensitivity obtainable 
with the proposed technique. We plan nevertheless to perform 
detailed modeling of atmospheric effects using physical op- 
tics, as has recently been done with reference to LGS ll2i . 



VI. CONCLUSIONS 

In conclusion, we have presented a promising alternative to 
satellite missions for the measurement of geomagnetic fields. 
The proposed method requires only ground-based apparatus, 
and is consequently substantially less expensive per sensor 
than a satellite formation, while still achieving high magne- 
tometric sensitivity. We anticipate that the low cost of de- 
ployment will make possible large-scale magnetic mapping 
and monitoring applications at the 100-km length scale, with 
temporal and spatial coverage that would be difficult to ob- 
tain by current techniques. Furthermore, as satellites can- 
not be operated as low as 100 km in altitude without exces- 
sive drag and heating, remotely-detected mesospheric mag- 
netometry promises superior spatial resolution of terrestrial 
sources. In addition, the technique offers to supplement exist- 
ing ground-based magnetic observatory data, allowing high- 
precision magnetic monitoring from a mobile platform with- 
out the requirement of a large-area, remote, and magneti- 
cally clean observation site on the earth's surface. We are 
currently constructing a 20-W-class laser projection system 
and working to implement this technique in proof-of -principle 
magnetic -field measurements. 

The authors acknowledge stimulating discussions with Pe- 
ter Milonni, William Happer, Michael Purucker, and Stuart 
Bale. This work is supported by the NGA NURI program. 
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